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ABSTRACT 



Aims. We study the photometric and structural properties of the star cluster system in the late type Sc spiral NGC 3370. 
Methods. BVI observations from the Advanced Camera for Surveys on board of HST are used to analyse in detail the colours, 
magnitudes and spatial properties of cluster candidates. The final catalogue of sources used for the study is composed by 277 objects. 
Results. The colour distributions of cluster candidates appear multi-modal. Although firm age constraints need the use of more age- 
sensitive indicators, the comparison of cluster candidate colours with the colours of Galactic and Magellanic Clouds star clusters, 
suggests an age difference between the various sub-peaks, with a red old sub-system, a rich population of intermediate age (~1 Gyr), 
and a blue tail of very young (below ~100 Myr) clusters. The luminosity functions appear normal for this type of galaxy, as for the 
distribution of cluster effective radii (Reff ). Our analysis suggests the presence of a peak in the R^ff distribution at ~3.0 pc, with 
blue (likely young) cluster candidates showing smaller radii respect to red (likely old) objects. Finally, inspecting the properties of 

candidate globular clusters, we find a colour distribution matching with Galactic Globulars, with a median [Fe/H] 1.5 dex, though 

a non negligible tail towards lower metallicities is also present. 

Key words, galaxies: individual: NGC 3370 ~ galaxies: star clusters 



1. Introduction 

The study of star clusters in galaxies provides an accurate and 
relatively straightforward tool to unveil the mechanisms that pro- 
duced the present organisation and evolutionary properties of 
stars in the galaxy. Star clusters allow to inspect the host galaxy 
both as a system and as single star clumps hosting -relatively- 
simple stellar populations. Moreo ver, stellar clusters a re the cra- 
dle of a large fraction of stars (e.g. lLada & Ladal2003L and refer- 
ences therein). As a consequence, understanding the properties 
of the Star Cluster (SC hereafter) system in a galaxy helps un- 
derstanding the properties of the host galaxy. Furthermore, SCs 
are relatively simple stellar systems - in first approximation sin- 
gle age and single metallicity stellar populations - much simpler 
than the stellar system of the host galaxies, providing the natural 
bridge from local single star analysis to distant unresolved stellar 
population studies. 

Motivated by the chance to analyse the evolutionary prop- 
erties of the host galaxy, in the last decade, especially af- 
ter the HST refurbishment, a great deal of work has been 
dedicated to the study the properties (luminosity, radius, 
mass, et c.) of th e SC sys tem in spir a l and irregular galax- 
ies (e.g. Whitmore & Schweizer 1995; Whitmore et al. 1999; 
Larsen & Richder ,2000; Bik et al.. .2003.; de Griis et al., ,2003^ 



Gieles et al ■I l2006allb[ just to quote few examples). 

In this work we present a study of the SC system in the spiral 
galaxy NGC 3370, based on archival multi-band images taken 
with the HST Advanced Camera for Surveys. 

^The galaxy is a nearly face on Sc spiral (Hubble type SA(s)c, 

Ide Vaucouleurs et al. 1991), showing an intricate tangle of arms, 
and mass similar to our own Milky Way . NGC 3370 was th e 
host of the type la Supernova SN 1994ae dvanPvket al.ll 19941) . 



Th e light-curve of th is SN la is one of the nine templates used 
bv lRiess et al.l (Il996h to introduce the SNe la multicolour light- 
curve shapes method to estimate distances. With the advent of 
the high efficiency ACS camera on board of HST, and thanks to 
the specific properties of the galaxy (face on, relatively nearby 
spiral), NGC 3370 has been targeted by an HST program to 
study the light-curves o f Cepheids, to pr ovide a direct calibration 
of the SN la distance (Ries s et ani2005l) . The program, carried 
out on optical V and 1 images, has been integrated by comple- 
mentary B band exposures for the Hubble Heritage program. 

The optical BVI data available offer an excellent opportunity 
to investigate accurately not only the photometric, but also the 
structural properties of the SC systems in this spiral, thanks to 
the high resolution and small PSF FWHM of the ACS. 

The paper is organised as follows. In section 2 we present the 
observational dataset used, the data analysis carried out, and the 
detection and selection criteria chosen for SC candidates. The 
study of photometric and spatial properties and the comparison 
with models is discussed in Section 3. We present our conclu- 
sions in Section 4. 

2. Observations, Data Analysis and Star Cluster 
selections 

The ACS/WFC camera observations of NGC 3370 used in this 
work come from two different proposals. The proposal ID #935 1 
(RI. A. Riess) consists of very deep F555W (~ V) and F814W 
band (~ I) observations to study the Period-Lum inosity relation 
of the Cepheids in the galaxy ( Riess et ani2005 ), with the pur- 
pose of deriving the distance of NGC 3370 and provide a precise 
luminosity calibration of the Type la supe rnova observed in the 
galaxy. The study realised by iRiess et al] provided the farthest 
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Fig. 1. ACS B, V and I composite image of NGC 3370, the image size is roughly 3.5 x 2.5 arcminutes (Credits: NASA, The Hubble 
Heritage Team and A. Riess - STScI). 



Table 1. NGC 3370 properties and observations 



RA & Dec (J2000) 
Morphological Classification 
Morphological T-type 
Total corrected B magnitude 
E(B-V) 

Distance Modulus 
F435W Exposure time (s) 
F555W Exposure time (s) 
F814W Exposure time (s) 



0h47m04.0s +17dl6m25s 
Sc" 

5.0± 1.1" 
~ 11.7" 
0.031'' 

32.29 + 0.06'^ 
9600 
57600 
24000 



" Data from Hyp erleda (http://leda.univ -lyonl.fr/) 
^ Extinction fromlSchlegel etal.Ul998l) 
' iRiess etal.l ( [200a> 



direct measurement of Cepheids. The other set of F435W band 
(~ B) images was obtained for the Hubble Heritage (proposal ID 
#9696, P.I. K. Noll). A colour image of the galaxy is shown in 
Figure [1] The main properties of NGC 3370, and of the obser- 
vational dataset us ed are revie wed in Table [T] For the distance 
modulus, we adopt lRiess et al.l value, /jq - 32.29 + 0.06, derived 
from an empirical calibration of Cepheids' Period-Luminosity 
relation, and an LMC distance modulus of 18.50 mag. 

Standard reduction was carried out on raw images usin g 
the APSIS image processing software (iBlakeslee et al.l l2003l) . 
Since ACS/WFC images are under-sampled for wavelengths at 
A $ 11,000 A, and thanks to the fine dithering of the input 
images, a resolution better than the ACS detector pixels size, 
0.05"/pixel, can be achieved with these images. We adopted a 
resampling of 0.035"/pixel as input parameter for APSIS, which 
corresponds to ~ 5 pc/pixel at the galaxy distance. 

After cosmic-ray rejection, alignment, drizzling, and final 
image combination, we derived the sky background from the 
median counts in the image corner with the lowest number of 
counts. This assumption is supported by the fact that the sur- 
face brightness profile of the galaxy is practically constant for 
r> 120" (Fig.H. 

The source detecti on was obtained using SExtractor 
jBertin & Arnouts !'l996) independently on the V and I frames. 
Due to the low S/N of B-band images, the source detection in 
this band was done using SExtractor in the pixel association 
mode (ASSOC mode), adopting as reference the coordinates of 
sources detected in the I-band frame. We used a minimum of 5 
connected pixels and 2cr threshold for sources detection. This 
resulted in a list of ~20,000 objects. The coordinates of sources 
detected from V and I frames were cross-checked, allowing a 
maximum wandering of 4 pixelfl then we performed a first 
very bland selection of SC candidates based on the SExtractor 
parameters star/galaxy classification, object major-to-minor axis 
ratio, and FWHM (object with CLASS.STAR = 0.00, elon- 
gation AJMAGE/BJMAGE > 2.5 and FWHM>12 were re- 



' The alignment between single HST frames is generally better than 
the maximum wandering adopted here. Our choice has been motivated 
by the large spatial variability of NGC 3370 surface brightness profile, 
and its wavelength dependence, which may lead to dilferent estimates of 
the source centroid in different bands. A posteriori, we verified that for 
none of the SC candidates in the final sample the coordinate wandering 
between V and I frames exceeds ~2 pixels. 




50 100 150 200 

r = Vab (arcsec) 

Fig. 2. Upper to lower panel: radial I, V and B band surface 
brightness profiles of NGC 3370. The linear fit to the data at r $ 
100"is shown with a dotted line. The horizontal dotted line at 
r>120"shows the background level in the three observed bands. 



jected). Moreover, we rejected the sources fainter than m/=27.5, 
and mv-28.5 ma^ Using such limits we ended up with a list 
of ~8000 SC candidates whose structural properties were then 
analysed. 

Given the linear pixel scale of these NGC 3370 images, and 
the ACS Point Spread Function FWHM of ~2.5 pixels, we are 
in the condition of using the i shape software (Larsen 1999|) 
to estimate the spatial extension of SC candidates. This soft- 
ware package allows to determine the physical characteristics 
of SC candidates {Rejf , elongation, etc.) by convolving analytic 
profiles with the PSF, then best-fit parameters are obtained by 
the minimisation. Ishape provides accurate intrinsic FWHM 
values of slightly resolved objects down to 1/10 the FWHM of 
the PSF, corresponding to ~1.2 pc in our case. In other words, 
provided that a good PSF sampling can be obtained, effective 
radii of star clusters ^ 1.2 pc can be retrieved. The input PSF 
to be used for ishape needs to be subsampled by a factor 10 
with respect to the input science images. As suggested by the 
i shape handbook, we generated our PSF using the SEEPSF task 



' The aim of the selection based on the SExtractor morphological pa- 
rameters CLASS_STAR, and AJMAGE/BJMAGE is to remove from 
the candidate list all obvious non SCs. Visual inspection of objects re- 
jected in this way revealed that they are mostly extended, irregular, dif- 
fuse or elongated objects. Moreover, in a further selection based on 
ishape output parameters (see end of this section), we rejected candi- 
dates with elongation estimated by ishape >2. The limiting magnitudes 
adopted here are, instead, related to the S/N ratio of ishape. Tests made 
on sources fainter than the limiting magnitudes adopted showed that the 
S/N is in all cases much low er than the minimum acceptable value ~30 
suggested bv lLarsenI ( Il999h . 
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Fig. 4. Some examples of visual classifications codes. For each 
source we show the original I-band image (upper sub-panels) 
and the ishape residuals (lower sub-panels). Objects associated 
to the classes 33 (large residuals, crowded area), 100 (possible 
background galaxies), and stars are rejected in the final catalogue 
of SC candidates. 



Fig. 3. Central panel: effective radii from the I-band frame versus 
the V-band R^ff (dots). The full squares show cluster candidates 
whose Rgff agree within 2cr in both bands (see text for details). 
The lower and left panels show the histograms of radii for the 
whole sample of sources (dotted lines), and for the SC candi- 
dates with matching radii (solid line). The dashed lines show the 
1 .5 pc and 20 pc cut radii. 



of DAOPHOT within IRAF0. To avoid contamination from clus- 
ters, the list of stars to determine the PSF was selected on the 
basis of the FWHM, colours and magnitudes from DAOPHOT. 
We finally used 10 bright stars (all common to the catalogue of 
comparison stars in Table 3 of Riess et al.), evenly distributed 
on NGC 3370, and all at distances > 50" from the centre of the 
galaxy. It is useful to anticipate here that a posteriori all the ob- 
jec ts used as te mplate for the PSF, and all the Cepheids detected 
bv lRiess et alj were catalogued as stars by ishape. 

Among the various profiles available with ishape, we 
choose a Moffat profile with a power law index -2.5, since it 
gave generally better residuals compared to other choices. A 
Moffat profile has already been shown to be effective in fitting 
the profiles of both young and old star clusters, though with 
slightly smaller exponents with respect to our c hoice (see, e.g., 
lElson et alJll"987l: lLarsenl|l999l: iMackev & Gilmor e 2003). As 
discussed in iLarsenld 19991) . the details of the model profile cho- 
sen are not an issue as long as the the Re/f values are derived 
for sources with intrinsic sizes smaller or similar to that of the 
PSF, as in the present case. The analysis was carried out inde- 
pendently on the I, V, and B frames, although only I and V band 
data were finally used to estimate the average effective radii of 
SC candidates. 

The ishape output parameters were used to further con- 
strain the list of SC candidates. In particular, we reject the ob- 
jects: 



IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



- whose effective radii Re/f differ more than 2cr in the two 
bands (Figure [3]l. The cr are derived using standard error 
propagation on ishape fitted parameters and uncertainties. 
In few cases the output of the fit is very good in both bands, 
so that the rms is small for V and I Re/j- , and the object will 
end up rejected by a simple 2cr rejection criterion. However, 
if the systematic uncertainties arising, for example, from the 
PSF, the analytic profile chosen, etc., are taken into account, 
then the two sources should be included in the catalogue. 
After several tests, we decided to use a default systematic 
uncertainty of ~15% Re/f ■ Then, for the 2cr selection we 
used the maximum cr between the rms from ishape, and 
-15% R,ff ; 

- with S/N,j,„p, < 30 (lLarsenllT999l) : 

- with too small/large Re/f . The specific properties of Moffat 
profile chosen allow to estimate effective SC radii be- 
low the ~ 1 .2 pc limit aforementioned (see iLarsenl 1 19991 
and ishape user's guide). However, we choose to adopt as 
lower and upper rejection limits Rg/f < 1 .5 pc, and Rgjf > 
20 pc, respectively. As a re f erence , using the data from 
iMcLaughlin & van der Marell (l2005l) . ~ 1% of massive star 
clusters in the Milky Way and its satellites have R,,// <1.5 
pc, while ~ 3% have Rgjf > 20 pc; 

- elongation (major to minor axis ratio) > 2.0, comparable to 
the maximum observed in Mag ellanic Clouds star clusters 
dvan den Bergh & MorbevlI984 . 

With these selection criteria we obtained a list of ~630 SC 
candidates. It is worth to emphasise the strongest constraints 
come from the S/N limit, and the Re/j- (most of the ~ 8000 
sources are actually classified as stars, see Figure[3]l. 

To verify the goodness of the ishape best fit parame- 
ters we have visually inspected all sources, introducing a sub- 
classification of sources depending on the residuals of SC candi- 
dates subtraction, and on the crowding. The visual classification 
code was used for further selections. Figure |4] shows some ex- 
amples of sources with different visual classification. The visual 
class code that we used is a number composed by two digits, the 
first associated to the goodness of ishape residuals, the second 
to the level of crowding. For both digits the range adopted is 1- 
3. For the first digit 1/2/3 means good/satisfactory /large residu- 
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Fig. 5. Aperture corrections in B, V, and I-band versus the SC 
candidate radius. The solid lines show the best fit linear equation. 



als, respectively. The second is related to crowding with 1/2/3 
meaning no crowding/less/more than three sources within 20 
pixels from the SC candidate. As an example visual class - 
12(31) means good residuals and 1 or 2 sources within 20 
pixels (large residuals and no other source within 20 pixels). 
A default visual class - 100 is adopted for diffuse objects 
(likely background galaxies) showing very large residuals in the 
i shape subtracted frame (see Figure|4|i. 

The aperture photometry of all the sources selected has 
been obtained using the Phot/Apphot task under IRAF. Due to 
the high, strongly variable background in some regions of the 
galaxy, we made different tests using 3, 5 and 7 pixels aperture 
radius, and 10, 16 pixel annulus for background determination. 
As a final choice, we adopted 3 pixels aperture for the photome- 
try, while the background was estimated in an annulus 16 pixels 
inner radius and 3 pixel width. 

Dealing with slightly extended objects having different radii, 
the aperture coiTection depends on the radius of the cluster As 
a consequence, we fitted the aperture correction versus R,,// . To 
avoid contamination from galaxy background, uncertain Rgjf es- 
timates, and crowding, we decided to use only sources brighter 
than mi^ipix - 26.5 mag, i shape S/N>50, visual class < 23, and 
local surface brightness background jiy > 22.5 mag (galactocen- 
tric radius Rgc ~ 35"). Aperture magnitude was derived in 3 and 
6 pixels (~ 15 and 30 pc) apertures over the selected sources. 
Figure |5] shows the aperture correction derived as the difference 
between 3 and 6 pixels aperture magnitudes, versus the effec- 
tive radii. The best fit linear correlations shown in the figure for 
each passband was then used to coiTect the magnitudes of all SC 
candidates in B-, V-, and I-band. 

The correction from 6 pixels to "infinite" radius, as well as 
the transformations from the F435W/F555W/F814W ACS filter 
system to standard B/V/I filters, and f oreground extinction, were 
applied following the prescriptions bv lSirianni et alj (l2005h . 

Once all corrections were made, we added new constraints 
to further clean the final list of SC candidates. In particular, to 
eliminate objects with very high internal extinction and possible 
distant background galaxies, we considered only those sources 
with V - I < 1.5 mag, B - V < 1.2 mag (typical maximum 
colours for globular clusters), and a maximum photometric un- 
certainty A{V - I) < 0.15 mag. SC candidates showing large 
residuals and in crowded regions, candidate background galax- 
ies (visual class =33 and 100, respectively) as well as sources 
located in regions with strongly variable background (to avoid 
objects with photometry significantly contaminated by internal 
extinction) have been rejected, too. The final list of SC candi- 
dates matching with the required constraints is composed by 277 
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Fig. 6. Colour-colour diagram of the SC candidates. Dark-grey 
dots, empty and full circles mark the positions of the sample 
after the first, second and third (final catalogue) selections de- 
scribed in Table |2] Th e grey shaded area marks the region of 
SPoT-SSP models by Raimondo et al.l (2005). The aiTow shows 
the direction of the reddening vector Five pointed stars are 
"outlier" SC candidates, possibly associated to intermediate age 
clusters. The insert panel shows the SC candidates (full dots), 
the SC outliers (five-pointed stars), and the standard SPoT SSP 
models (grey-shaded area). The dark-grey area shows the re- 
gion occupied by S SP realisations of the recent SPoT models 
by Raimondo (2009), computed paying particular attention to 
the stochastic effects of TP-AGB stars. Only ages between 0.3-2 
Gyr and [Fe/H] 0.4dex are plotted (see text for details). 



objectfl A schematic view on the effects of the various selection 
criteria applied, with more details on the numbers, is reported 
in Table |2l The table gives the number of objects before (A^,„) 
and after (Nout) each selection criteria is applied, together with 
the list of criteria adopted and the number of candidates rejected 
with the chosen parameter (Nrej)- It is worth noting that the same 
object might satisfy two or more rejections, thus the total number 
of objects rejected, Nrepot, is not simply obtained by summing 
up the numbers in the Nrej lines. The positional, photometric and 
structural properties of the final catalogue of sources are reported 
in Tables 

3. Discussion 

3.1. Analysis of Colours 

Figure |6] shows the colour-colour diagram of the final catalogue 
of SC candidates. In the figure it is shown also the effect of the 
various steps of the selection. Dark-grey dots, empty and filled 
circles mark the position of the sample after the selections based 
on SExtractor, i shape , and photometric parameters, respec- 
tively. The full circles in the figure refer to the final selected 
sample of SC candidates. Also shown in the figure, with a grey 



The complete catalogues of sources can be retrieved upon request 
to the authors. 
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Table 2. Number of objects rejected per selection criteria 
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First selection 
Nin = 17767, 
Sel. Criteria 


: SExtractor selection criteria 
A^„„, = 8081, A',,,,„, = 9686 
CLASS_STAR=0.00 Elong.^^,, >2.5 
675 1652 


nil > 27.5 
3251 


niv > 28.5 
3962 


FWHM, > 12 FWHMv > 12 
2452 3287 


Second select 
A';,, = 8081, A 
Sel. Criteria 

A' ■ 


ion: Ishape selection criteria 
'„„, = 635, N„„„, = 7446 

2aR,ff Elong.,s,,. > 2.0 
1414 1293 


S/N<30 
5051 


Reff > 1.5 pc 
5277 


Reff < 20 pc 
953 


Third selectio 
Nin = 635, A^,, 
Sel. Criteria 

A' ■ 


n: photometric selection criteria 
„, = 277, N„jj„ = 358 

A(V-/)>0.15 vis. class>33 
17 286 


V-I> 1.5 
66 


B-V > 1.2 
19 





Table 3. Structural and photometric properties of selected SC candidates. [See the electroniv version of the Journal for the complete 
version of the Table.] 



ID RA(deg) Dec(deg) Rg,(") R,ff (pc) ± my ± B-V ± V-I ± v.c. 



5 


161.7465 


17.2766 


71.65 


2.04 


0.25 


25.39 


0.01 


0.23 


0.02 


0.28 


0.01 


21 


10 


161.7527 


17.2907 


79.03 


5.50 


0.45 


26.15 


0.02 


0.67 


0.06 


0.94 


0.03 


12 


13 


161.7513 


17.2858 


69.83 


4.68 


0.92 


25.40 


0.02 


-0.04 


0.03 


-0.09 


0.05 


23 


14 


161.7552 


17.2964 


91.56 


5.58 


0.41 


25.98 


0.01 


-0.03 


0.03 


0.00 


0.03 


22 


15 


161.7556 


17.2980 


96.25 


13.67 


0.21 


24.57 


0.01 


0.51 


0.01 


1.30 


0.01 


21 


17 


161.7537 


17.2915 


79.24 


6.86 


0.45 


25.77 


0.01 


0.58 


0.04 


0.82 


0.02 


23 


18 


161.7527 


17.2879 


71.40 


6.62 


0.18 


24.46 


0.01 


0.28 


0.01 


0.58 


0.01 


22 


21 


161.7506 


17.2792 


60.20 


6.35 


0.44 


26.26 


0.02 


0.22 


0.04 


0.69 


0.03 


12 


22 


161.7523 


17.2850 


65.42 


2.13 


0.38 


25.09 


0.01 


-0.09 


0.02 


0.03 


0.01 


22 


23 


161.7545 


17.2898 


72.80 


4.50 


0.34 


25.69 


0.02 


0.70 


0.06 


0.86 


0.03 


23 


24 


161.7519 


17.2828 


61.81 


4.00 


0.45 


25.76 


0.02 


0.17 


0.04 


0.56 


0.02 


23 


25 


161.7530 


17.2856 


64.93 


3.84 


0.40 


25.58 


0.01 


0.38 


0.03 


0.56 


0.02 


12 


29 


161.7495 


17.2750 


60.90 


3.89 


0.46 


26.04 


0.01 


0.30 


0.03 


0.40 


0.02 


22 


30 


161.7550 


17.2895 


70.74 


13.23 


0.46 


26.17 


0.02 


0.34 


0.06 


1.38 


0.03 


12 


33 


161.7529 


17.2837 


60.83 


3.06 


0.10 


25.40 


0.01 


0.18 


0.03 


0.39 


0.02 


23 


34 


161.7573 


17.2952 


84.98 


2.54 


0.08 


24.64 


0.01 


0.70 


0.02 


1.05 


0.01 


22 


36 


161.7603 


17.3021 


105.49 


2.74 


0.57 


25.90 


0.01 


0.42 


0.02 


0.39 


0.01 


22 


37 


161.7549 


17.2870 


64.02 


5.79 


0.24 


26.37 


0.02 


0.62 


0.07 


1.16 


0.03 


13 



area, is the region of the SPoT Simple Stellar Populations mod- 
els for metallicity between [Fe/H] ~ -2.5 0.3 dex, and a ges 
in the interval 25 Myr -H 14 Gyr (see iRaimondo et ani2005l for 
more details on SPoT SSP models). As expected, the data over- 
lap to SSP models, since SC are, at least in first approximation, 
well represented by SSPs. It must be emphasised, however, that 
such overlap is not a trivial result. In fact, in the selections de- 
scribed in the previous section, the only constraint on colour was 
on the maximum V-I and B-V allowed. 

One interesting feature appearing in Figure |6] is the non- 
negligible number of objects located at B-V ~0.4 and B-/ ~1.6 
(five pointed stars in the figure). After different counter-checks, 
we concluded that there is no obvious reason to reject these clus- 
ters from the final catalogue. As shown by the arrow in the pan- 
els, even an internal extinction E(B-V)~ 0.5 mag would not shift 
these SC candidates in the region of SSP models. Moreover, such 
high extinction is unlikely since such objects do not appear ob- 
viously associated to regions of high internal extinction, while 
they are rather evenly distributed over the area of the galaxy. 
One possible interpretation of these objects is that they are in- 
termediate age SCs (0.3 < t{Gyr) < 2). Within this age inter- 
val, the emission of Thermally Pulsating AGB stars (TP-AGBs) 
dominates the light emitted by the SC, in particular in the long- 
wavelength regime ( ^ 8000A; see, e.g. Maraston 2005, and ref- 
erences therein). However, the presence of stars in such evolu- 
tionary stage is strongly influenced by stochastic effects. As a 



consequence, the stochasticity due to the number of cool TP- 
AGBs can significantly affect the integrated colours of the stellar 
system. To highlight such behaviour, the insert panel in Figure|6] 
shows the sample of selected SC candidates, and the "outliers" at 
B-V ~QA and B-/ ~ 1.6 mag (five-pointe d stars), together w ith 
the region where the recent SSP models bv lRaimondd(l2009l) for 
intermediate ages and [Fe/H] ~ -0.4 dex are distributed (dark- 
grey area). The new set of models used has been computed with 
an updated version of the SPoT code, paying particular attention 
to the statistical effects generated by Thermally Pulsating AGBs. 
In the figure, the dark-grey area shows the distribution of mod- 
els when the stochastic effects due to TP- AGB stars are properly 
taken into account. As shown in the insert, these models predict 
the existence of the outliers (optical models have been provided 
as a private communication by G. Raimondo, for more details on 
models and near-IR colours see lRaimon dol2009h . In a following 
section we will discuss more in details the data to models com- 
parison for the whole sample of SC candidates. 

Figure [T] shows the colour magnitude diagrams and colour 
histograms for the SC candidates. It is not surprising that the 
whole population of SC candidates in NGC 3370 does not have a 
unimodal colour distribution, showing the presence of sub-peaks 
possibly associated to different SC sub-populations. 

As a matter of fact, since NGC 3370 is a late type Sc spi- 
ral, it is expected that this class of galaxies is actively forming 
stars, and star clusters. Thus, populations of young/intermediate 
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Fig. 7. Left panels: B-V colour magnitude diagram (upper panel) 
and colour histograms (lower panel) of SC candidates. Right 
panels: as left, but for B-I. The arrows indicate the median colour 
ofGGCs. 




Fig. 8. Colour histograms of the comparison galaxies (as labelled 
in left panels). Dotted lines show NGC 3370 SC candidates. For 
NGC 3367 and SMC a different binwidth is used, because of the 
smaller number of SCs (see Nobj labels in each panel). 



age SC are expected to host in the galaxy, in addition to the old 
population of Globular Clusters (GCs). To highlight this point, 
in the colour histograms shown in Figure |7] we have added 
an arrow showing the median colours of Galactic GCs (GGCs, 
(B-y)o = 0.73 +0.12, (B-/)o = 1.65 +0.12 mag for -90 GGCs 
with My ~ - 5.0 mag, data from Harris 1996). It is rather ob- 
vious, then, to conclude that the red peak observed in the lower 
panels is actually due to the GC population in the galaxy. 



If the intermediate peak and the blue tail of SC candidates 
are taken into account, still using as reference the system of 
young open clusters (OCs) in the Galaxy, we find that the OCs 
at B - y -0.35, and B-I -0.85 mag have ages log t(yr) -8.4 
(based on only 3 Galactic OCs), while blue OCs at B - V ~- 
0.15, and B-I — 0.25 mag have on average logf(3'r) ~ 7.0 
(based on 1 1 OCsfl Given the results of such comparison with 
Galactic star clusters, it is reasonable to conclude that the colour 
peaks observed in Figure|2]are mainly driven by age differences 
between the blue/intermediate/red SC candidates. 

The situation, however, might be complicated by the possi- 
ble presence of residual internal extinction. Since the galaxy is 
oriented nearly face on, the effect of internal extinction is min- 
imised. Moreover, as discussed in previous section, to reduce 
possible extinction present in the dusty, star forming regions we 
rejected the sources located in areas of strongly varying back- 
ground, or where dusty patches can be recognised. The pres- 
ence of the colour peaks discussed above is an additional indi- 
rect probe against a possible strong contribute of internal red- 
dening: a significant amount of dust extinction, in fact, should 
smear out such features. However, it is not possible to exclude 
at all the presence of some residual internal extinction. To have 
other indications about the possible nature of the colour range 
and sub-peaks in NGC 3370, we compare the B-I and B-V 
colour histograms of SC candidates with the same data for other 
galaxies, some of which are corrected for internal redden ing. In 
particular, we use the dat a from 1 ) Milky Way's GCs dHarrisI 
""6), and Open Clus ters jLata et alj |2"002l): 2) NGC 3367 from 
:fa-Barret oetan ( l2007h : 3) the data of ~ 20 S pirals from 
Larsen ( 1999); 4) M 33 data corrected for reddening ( Park & Lej 
2007); 5) SMC St ar Clusters from Rafelski & Za ritsky (2005|); 
6) LMC data from lBica et alj (Il996h . All the data used are in (or 
are transformed to) the standard photometric BVI bands. The 
comparisons between different galaxies are shown in Figure [8] 
and|9l In order to compare only SCs in the luminosity range of 
the SC candidates considered here, we have excluded all SCs 
having total absolute magnitude My ~ - 5 mag. 

Figure [8] shows B-I and B-V colour histograms of the 
SC in the MW, NGC 3367, the spirals from Larsen (1999), M 
33 and the SMC. The galaxies are ordered according to their 
morphological T-Type, all histograms are normalised to have 1 
in the most populated bin. In each panel of the figure the colour 
histograms of SC candidates in NGC 3370 are shown with dotted 
lines. 

By inspecting the colour histograms it is easy to reveal sim- 
ilarities and differences between the system of SC candidates in 
NGC 3370 and in other galaxies. One example is represented 
by the red sub-peak, which appears quite similar to the red sub- 
peak in the colour distribution of SMC clusters. Using the age 
estimates by iRafelski & Zaritskvl (l2005h we find that the ~ 20 
SMC clusters with B - / ^ 1.6 mag have a median age f ~ 5 
Gyr The range of colours generally shown by other systems is 
quite similar to that of NGC 3370. The only exc eption is repre - 
sented by the colour histograms of the SCsfrpmlLirsi^ (Il999l) . 
However, the selection criteria used b v iLarsenl exclude SCs with 
B -V > 0.45, which represents ~45% of our sample of objects. 
In addition, the very blue/young stellar clusters are rejected by 
Larsen, adopting a lower limit to the Ha - R colour Thus, we 
are lead to conclude that, ra ther than a systematic difference be- 
tween our data and iLarsenf s. the observed differences are due 
to the intrinsically different SCs selection. In any case, the peak 



' Ages and colours for OCs are from the WEB DA database, available 
at |http://www.univie.ac.at/webda, 
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Fig. 9. B-V colour histograms of LMC Star Clusters (solid lines) 
and NGC 3370 (dotted lines). The vertical blocks refer to LMC 
clusters selected according to the SWB class: the first, second 
and third row refer to LMC clusters with SWB> 0, 3, 6, respec- 
tively (see labels in upper panels). Horizontal blocks assume dif- 
ferent average internal extinctions for NGC 3370 SC candidates: 
first, second and third row have {E(B - V)i„,enmi) - 0.0, 0.1, and 
0.2 mag, respectively (see label in leftmost panels). 
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Fig. 11. Annuli division used to compute completeness func- 
tions. The regions are characterised by their median surface 
brightness: ju*' <20, juf between 20H-23, /if between 23^-26.5 
and fi*^ > 26. 5mag/ arc sec^. Each annulus is labelled and iden- 
tified with a different shade of grey. Also the loci of SC candi- 
dates are shown with five pointed stars. No SC candidate lies in 
the innermost annulus (#1, white area). 
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Fig. 10. Upper panels: colour histograms of SC candidates, the 
position of the three relative peaks is emphasised in each panel 
with dashed lines. Lower-left panel: spatial distribution SC can- 
didates selected according to their membership to the colour sub- 
peaks (see text). Five pointed stars/full squares/circles mark can- 
didates in the blue/intermediate/red tail respectively. 



of the distributions of'Larserfs data appears similar to the main 
peak of NGC 3370 SCs, though slightly shifted towards bluer 
colours. 

The SCs data of NGC 3367 - the galaxy of the sample mor- 
phologically more similar to NGC 3370 - do not show obvious 
multi-modal colour distributions. The range of colours, and the 
position of the main peak are roughly similar to NGC 3370. As 
for what concerns SMC clusters, the colour distributions span 
over the same range, and the peak of SC sub-systems seems to 
be recognised here as in NGC 3370. 

We provide in Figure|9]a separate comparison to LMC data. 
The reason for keeping LMC data aside is that, in this case, 
we can play w ith the age of the selected clusters by using 
the SWB class jSearle et al.lll980h of LMC clusters, and also 
make a test on the internal extinction in NGC 3370. The fig- 
ure shows the B-V colour histograms of LMC clusters, cor- 
rected for internal extinction, and NGC 3370 SC candidates. 
Each row of panels in the figure shows NGC 3370 cluster can- 
didates dereddened adopting a different average internal extinc- 
tion, {E{B - V)iniernai) - 0.0, 0.1 and 0.2 mag from upper to 
lowe r panel, resp ectively. In addition, assuming the SWB classes 
from iBica et aO [l996,) . each column of panels is obtained as- 
suming a lower limit to the SWB. Specifically, the first column 
includes all SWB classes, while the other two include only LMC 
clusters with SWB>3 and 6, respectively (labels in the upper- 
most panels). All panels show NGC 3370 colour histograms 
(dotted lines) shifted according to the internal reddening as- 
sumed, and the LMC clusters colour histograms using the SWB 
selection criteria (solid lines). Some hints on age distributions 
and, possibly, internal extinction of the SC candidates can be 
derived by a careful inspection of this figure. The upper panels 
show that if old LMC clusters are taken into account (SWB> 6, 
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t ^ lO' yr), there is a good match between NGC 3370 red peak 
and the old clusters in LMC. Thus, we have an additional indica- 
tion that this population is, in fact, composed by old GC candi- 
dates. It is worth to notice, in the case of {E(B - V)intenmi) ~ 0.1 
and SWB> 3 (central panel in the figure), the similarity be- 
tween the main colour peaks of the two galaxies. Clearly, be- 
cause of the differences between the two galaxies, it is not pos- 
sible to conclude whether or not such good matching is a con- 
sequence of an average internal extinction of ~0. 1 mag on a SC 
system in NGC 3370 with an age distribution similar to LMC. 
It is, however, reasonable to take this value as an upper limit 
for {E{B - V)in,emai)- Furthermore, the lower panels, which as- 
sume larger internal extinction, show that such average redden- 
ing would move the main peak to a colour corresponding to ages 
of the order of logf(yr) ~ 7.5 - obtained taking MW OCs and 
LMC clusters ages as reference. Although young SCs are ex- 
pected to be present in late type, actively star forming galaxies, it 
is not reasonable to conclude such predominance of very young 
objects. The position of the red peak is not an issue in the latter 
case, since it is expected that old clusters are, on average, less 
affected by dust extinction being far from actively star-forming 
regions. Thus, the shift of the red peak for {E{B-V)i„,ernai) - 0.2 
mag is not probable. However, since this test is carried out us- 
ing a common average value for internal extinction, we cannot 
rule out high extinction values on single bluer/younger cluster 
candidates. 

In conclusion, the comparison of the colour histograms for 
NGC 3370 SC candidates with other SCs from literature shows 
that the range of colours appear normal with respect to simi- 
lar class of galaxies. Clearly, the comparisons shown cannot be 
used to constrain ages, or any residual internal reddening, be- 
cause of the intrinsic differences existing between the galaxies. 
However, we can give some limits, in particular, it is reasonable 
that NGC 3370 hosts three sub-populations of SC candidates: 
one with colours typical of old GCs, one with very blue colours, 
characterising a population of very young SCs, and a rich pop- 
ulation of intermediate age SC candidates. This expectation is 
partly supported by the spatial distribution of the SC candidates 
in the three colour peaks, shown in Figure [10] In the figure we 
used the position of the three relative maxims to define a prelim- 
inary centre of each sub-peak in B-I and B-V. Then, the median 
and standard deviation of the sole SCs within half-distance from 
the sub-peaks is obtained (dashed-lines in the upper panels). All 
SC candidates within 2cr from the median value of the peak in 
both colours are then selected. The SCs selected in this way are 
shown in the lower- left panel of the figure: five pointed stars/full 
squares/circles mark candidates in the blue/intermediate/red tail 
respectively. By inspecting the spatial distribution of clusters, we 
recognise that the SC candidates in the blue tail are mostly lo- 
cated near the galactic centre, while the star clusters in the other 
colour bins are more spread over the galaxy, with the ones in 
the red peak being on average more distant from the centre of 
the galaxy, as expected for GC candidates. The median effec- 
tive galactocentric distance of each SC class is estimated and 
reported in the figure. 

3.2. Luminosity Functions 

Starting from the work by lElson & Fall (Il985h on LMC SCs, 
empirical data from a large sample of spiral and irregular 
galaxies have shown that the Star Cluster Luminosity Function 
(SCLF) usually follows a power law, with a sharp turnover at 
faint magnitudes. In this section we analyse the LF profile of the 
SC candidate system in NGC 3370. 
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Fig. 12. B, V and I-band completeness functions for the each of 
the four annuli considered (labelled in the panels). For each one 
of the 4 annuli, we have carried out several numerical experi- 
ments shifting the grid of artificial stars. Different experiments 
are shown with different line types in each panel. The horizon- 
tal dotted line marks the 90% completeness limit in the annulus. 
The rightmost panels show the median completeness function 
for each annulus and filter 

To properly study the SCLF, we need to take into account the 
completeness corrections in the different bands. However, be- 
cause of the selection criteria adopted, this is a complex task. In 
our case, in fact, the completeness functions should take into ac- 
count /) the photometric completeness, //) the selection on sizes, 
and ///) the visual identification. While the first two contributions 
to the completeness can be somehow obtained from numerical 
experiments, the third one is less straightforward to derive, if fea- 
sible at all. Moreover, the completeness fraction is different for 
objects of different sizes. For such reasons we choose to adopt 
the following approach: the galaxy is divided in four concentric 
regions, selected on the basis of the median surface brightness 
of the galaxy. In particular, we choose fif < 20, 20 < juf < 23, 
23 < fi*^ < 26.5, and fi*'^ > 26.5 mag/arcsec^ as separation be- 
tween the regions #1, #2, #3, and #4. Figure [TT] shows the four 
regions and the position of SC candidates. In each of the four an- 
nuli, then, we run completeness tests separately. To add artificial 
SCs we used DAOPHOT/ADDSTAR, adopting a template SC 
candidate derived from the best candidates in our list. Accurate 
completeness functions would need different numerical experi- 
ments obtained with SC templates of various Re/f . We have not 
carried out such detailed analysis. It is worth to remark, though, 
that the median size of template SC, {Refj) ~ 3 pc, is represen- 
tative of the vast majority of the SC candidate population (see 
next section). 

To derive completeness functions, SC candidates at fixed 
magnitudes are artificially added in the frame on a grid of 100 
pixel separation both in x and y coordinates, the photometry of 
artificial sources is then obtained in the same way as for the 
original data. The comparison of the input and output photo- 
metric catalogues is used to evaluate the completeness fraction 
/ = Nsc deteciedlNsc added- The experiments are repeated for /) 
different magnitudes (between 15 and 30 mag, with steps of 0.1 
mag), and //) for different positions of the grid of stars. The re- 
sults of the completeness tests are shown in Figure[T2] The com- 
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Fig. 13. Left panel: B-band Luminosity Function of SC candidates. The dot-dashed line shows the power law fit to the data. The 
exponent of the fit is also reported in the panel. Middle panel: as upper panel, but for V-band data. Right panel: as upper panel, but 
for I-band data. 



pleteness fraction is quite low in the innermost annulus (#1), 
though we recall that no source in the final catalogue lies in such 
area of high, and highly variable galaxy background (FigurefTTli. 

Using the completeness functions so derived, we estimated 
the corrected LP of SC candidates, shown in Figure [13] The ob- 
served sharp falloff is an artifact due to the sample incomplete- 
ness at faint magnitudes. In the three panels of the figure we 
plot a fit to the data, the slope a obtained from a power law fit 
dN/dL oc L" is also reported. To minimise the effect of incom- 
pleteness, the upper limit of the fit is chosen so that the complete- 
ness fraction in the innermost useful annulus, #2, is / ^ 0.5 0. 
Notwithstanding the complex issue of completeness, the power 

law fit to the data provides exponents a 2.2 si milar to what 

typically found in spirals and starburst galaxies (e.g. Miller et 
[l997t IWhitmore et al.lll999t iLarsenlllOOl Ide Griis et al.ll2003l) . 
We also find, for all three passbands considered, steeper slopes 
on the bright sides of the LPs, with power law exponents as large 
as a ~ 3.0. In all cases, however, the new exponent agrees within 
Icr uncertainty with the single power- law slopes reported in the 
figure. The presence of a bend in the SCLF has alrea dy been 
found in various galaxies. As discussed in details by Giel es et all 
f2006aiih) . a SC system originating from a truncated mass func- 
tion (MF), possibly caused by a truncated MF in the giant molec- 
ular clouds, is expected to show such a bend, with steeper slopes 
on bright side of the LF. One other prediction of a truncated MF 
is the steepening of the LF at longer wavelengths. In spite of 
the uncertainties on the fitted slopes, this behaviour seems to 
be found in our data. Nevertheless, comparing maximum clus- 
ter luminosity (Mv,ma.K 9.5 mag) with the number of clusters 

brighter t han My ~ - 8.5 mag (N(Mv < -8.5) = 13), as pro- 
posed by iGieles et"an (2006a, see also the references therein), 
we find that the luminosity of the brighter cluster in NGC 3370 
can be accounted for by a size-of-sampl e effect, with a single 
power-law distribution (see Fig. 16 in Giel es et aDl2006al) . 

The consistency of the SCLF with a power law indicates 
no preferred scale for cluster luminosity (i.e. mass). However, 
"passive" evolutionary processes like tidal interactions, dynami- 
cal friction, and the "active" evolution of stars in the cluster are 
expected to alter the shape of the SCLF. On the contrary, the lu- 
minosity function of GCs (GCLF) in both elliptical and spiral 



galaxies shows a nearly univ ersal turn-over at My ~ -7.4 mag 
(iHarrisllIOOll; lRichtledl2003h . implying that present GCs have a 
preferred mass scale. By inspecting the data in Figure [13] we 
do not see any peculiarity near the regions of the GCLF peaks 
{my ~ 25 mag, and m/ ~ 24 mag), although this might well be 
a consequence of the poor number of GC candidates in our final 
catalogue (Sect. [33T l and the dominance of young objects. 

One more piece of information comes from the study of 
the correlation between the SCLF and Star Formation Rate 
(SFR) of the galaxy. As shown by different authors (e.g. 
Larsen & Richtler 2000; Billett et al. 2002; Larsen 2002), some 
properties of the SC system are correlated to the SFR, with 
the efficiency of cluster for mation being la rger in galaxies with 
higher SFR. In particular, iLarsenl (|2002|) analysed the maxi- 
mum luminosity expected from a random sampling of the LF 
(assuming a power law exponent a = -2.4) versus the SFR 
per unit area. The comparison of predictions to observations 
showed a satisfactory agreement. Using the LF shown in Figure 
[13] we find that the maximum observed SC candidate magni- 
-9.5 mag. Assuming the SFR=0.08 Mq yr 



tude is M'. 



by |Grosb0l & Dottorij 12008, a "crude calculation " based on the 
K-band properties of bright knots in the galaxy), and an area 
A~ 240 kpc^, we find that NGC 3370 data fall in the same obser- 
vational regions as the spirals in the Figure 16 of Larsen (200^, 
without any noticeable scatter with respect to other observational 
data. More over, the beha viour with respect to the predictions 
presented in lLarsenI (120021) . based on sampling statistics, is offset 
by nearly the same amount (~0.9 mag, or ~0.7 mag, depending 
on the model) as other spirals. 

3.3. About effective radii and spatial distribution 

In this section we analyse the Re/f versus various other charac- 
teristics of the SC candidates, or the host environment. Let us 
first inspect the Re/f distribution of the system. 

Upper panel of Figure [14] shows the distribution of R^ff us- 
ing both the radii obtained from the weighted average of V and I 
values (solid histogram), and those obtained using all BVI data 
(dotted line). The differences between the two distributions is 
negligible, mainly because the B-band R^jy have larger uncer- 
tainties, thus lower weights. As anticipated, we adopted the V 
and I averaged R^ff values. 



^ The average / between annuli #2, #3, and #4 at the faint limit is 
larger than 0.8 in all bands. 
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Fig. 14. Upper panel: R^ff distribution for SC candidates. Solid 
line histogram shows the mean Rgfj- obtained from I and V band 
data, the dashed-line shows the histogram obtained coupling B- 
, V- and I-band data. Lower panel: As upper panel but in log 
scales. The thick lines show a power law (only for Rgff ^ 4 pc, 
solid line), and a log-normal fit to the data (all R^jf , dotted line). 



The lower panel of Figure [14] shows the distribution of 
Reff in logarithmic scale. A power law fit to the data (solid 
line) for R^jf larger than ~3 pc, provides an exponent rj - 
-2.35 + 0.22. Such value is consistent with simila r results for 
other galaxies. For example, iBastian et al.l (l2005h found rj - 
-2.2 + 0.2 for the star clusters in M 51, although they fit with the 
same power law the whole set of Rgjf (2 pc < R^fj < 15 pc). In 
Figure [14] it is shown also that a log-normal distribution with a 
peak at ~3 pc (dotted line) provides a good fit to data. 

The median R^jj of the system is 4.2, with 3.0 pc standard 
deviation. If the three possible sub-populations of SC candidates 
(Sect. 3.1) are taken separately, we obtain as median and stan- 
dard deviation R.fj = 2.9 + 1.5, 3.8 + 1.8, 5.5 + 3.2 pc for 
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Fig. 15. Median R^fj estimated per colour bins. The slope of the 
linear fit to data and its nns are reported in the panels. 




Fig. 16. Average effective radius versus distance from the galaxy 
centre. 



the blue/intermediate/red peak, respectively (using only the SC 
candidates within Icr the colour intervals reported in Fig. [TO] 
the number of objects used are 7, 72, and 18 from blue to red 
peali[3)- Although the three median Rejf agree within Icr, we no- 
tice smaller R^jj for bluer SC candidates. Figure [T5]gives more 
details on this regard. In the figure we show the average Rgjf per 
colour bin for both B - I and Z? - V. In spite of the large scatter, 
it is possible to observe the presence of a non negligible correla- 
tion between the size and the colour of SC candidates, with red 
SCs showing on average larger Rejf respect to blue candidates. 
Such effect has already been observed in other SC systems (e.g. 
Scheepmaker et al. 2007). Again, if colour differences are asso- 
ciated to age differences, then the increase of Rgjf with colour 
points toward a dynamical evolution of clusters radii with age, 
with younger objects being more concentrated than older ones. 



We obtain R, 



■If 



3.3 ± 1.1, 3.8 ± 2.1, 5.3 ± 3.5 pc for the 



blue/intermediate/red peak using SC candidates within 2cr the colour 
intervals of Fig. 1101 
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Fig. 17. Average colour of SC candidates versus galactocentric 
distance. 



The spatial distribution of star clusters provides an additional 
piece of information on the properties and origins of the SC sys- 
tem and the host galaxy. One possible clue comes from the re- 
lation between the half-light radius and the galactocentric dis- 
tance Rgc, shown in Figure [16] Also in this case, the spread 
around the average R^ff for each Rgc bin is large, but a corre- 
lation over the entire Rgc range seems to be present. A linear fit 
to the data provides a s lope a ~ 0.30. By analysing the SC sys- 
tem of M 5 1 , Scheep maker et al.l (i2007) found a shallower slope, 
a = 0.12 + 0.02. However, similarly to M 51, we find no corre- 
lation, or a very weak one, between the colour of the SC can- 
didate and Rgc (Fig. [17). Finally, the spatial density of SC can- 
didates versus Rgc, shown in Figure [18] reveals that the surface 
brightness profile of the galaxy (in the figure with an arbitrary 
shift) agrees with the radial density profile of SC candidates. In 
the figure we omit the central ~2 kpc as SC candidates in such 
bright region are rejected. A linear fit to SC data provides a slope 
o^sc = -0.3 1+0.02, very similar to the ffsja/fl^//^/,, = -0.28+0.02 
from the surface brightness profile of the galaxy. Such behaviour , 
typically observed in late type galaxies (ISchweizer et al.|[T996l: 
iMiller et al.ll997l;lLarsenll999l) . might be an indication that field 
stars and star clusters have experienced similar dynamical relax- 
ation processes and, that star clusters at any location (and age?) 
are in dynamical equilibrium with the galaxy potential. 



3.4. Mass function from comparison with models 

Mass estimates of SCs are sensitive to the age of the cluster, be- 
cause the mass luminosity ratio changes with time, due to stel- 
lar evolution, and because of the dynamical evolution of clus- 
ter itself (tidal disrupti on, evapora t ion, e tc.). In the approxima- 
tion of no dissolution, Gieles et al ] (l2007h . using SMC data, find 
that the minimum detectable cluster mass scales with age f ac- 
cording to Mmin oc f^"^, or, in other words, the mass detec- 
tion limit increases by a fac tor of ~ 5 each age dex (see also 
iBoutloukos & Lamer sl l2003L for a similar analysis on clusters 
in the MW and nearby galaxies). As discussed in Section 3.1, 



O 




R,, (Kpc) 



Fig. 18. Radial density of SC candidates corrected for complete- 
ness (thick solid line). The completeness correction is only sig- 
nificant in the innermost two bins considered. The upper thin 
solid line shows the surface brightness profile of the galaxy 
shifted vertically by an arbitrary amount. The slope of a linear 
fit to both distributions is also reported in the figure. 



the present sample of clusters contains both massive old GCs, 
which have a peaked mass function, and young clusters, which 
form from a power-law mass function. Therefore, an age depen- 
dent mass-completeness affects our data. To properly constrain 
the mass of SC candidate, age sensitive data, like U-band pho- 
tometry, are essential, and a study of the SC masses drawn from 
only BVI data must be considered with extreme caution. 

Keeping in mind these warnings, in this section we tenta- 
tively derive rough mass estimates of SC candidates. To this pur- 
pose we use magnitudes and colours of SC candidates and com- 
pare them with stellar population synthesis models. In particu- 
lar, we take ad vantage of the specific ca pabilities of the Tera mo- 
SPoT models (ICantiello et al . 2003; Raimondo et al.ll2005lFl to 
obtain some hints on the mass distribution of the SCs. Drawing 
the mass distribution from a comparison of data to SSP models is 
a risky exercise, and could be misleading if one does not consider 
properly the limitations due to the assumptions done. The total 
mass of the S SP strongly d epends on the Initial Mass Function. 
We as sume d alScalol (Il998l) law; a diff'erent IMF, like a lSalpeted 
(Il955h or a lKroupal(l2001h IMF, would affect the total mass es- 
timates by few tens percent. However, that's rather unimportant 
for the rough, and relative analysis that we will carry out here. A 
second approximation comes from the fact that the age and mass 
associated to each SC candidate depends on the metallicity of 
the stellar cluster For this reason, we have assumed three metal- 
licities, representative of metal rich, intermediate, and and metal 
poor stellar populations, i.e. [Fe/H] ~-1.8, -0.7 and 0.0 dex, re- 
spectively. Then, we have analysed the three cases separately. 

To derive crude mass estimates of the SC candidates we used 
two different approaches, both based on the data to model com- 
parisons: first we compare SC candidates colour and magnitudes 



Visit the website www.oa-teramo.inaf.it/spot for details on models, 
code, and for downloads. 
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with SSP models for the three selected [Fe/H] separately, then a 
second estimate of masses is obtained in the colour-colour plane 
from median mass-luminosity ratios. 

For the specific purpose of this work we have computed 
SSP models assuming different total masses, i.e. 5.0 • lO"*, 1.5 ■ 
1Q\ 6.0■10^ 1.2-10\ 1.0-10<'Mo, and ages 0.03, 0.2, 0.5, 1.0, 
5.0, 10.0, 14.0 Gyr, for the three adopted metal contents. The 
upper panels in Figure [19] show a comparison of data to models 
for the selected [Fe/H] values, all computed ages and masses are 
shown. 

Assuming a defined [Fe/H] , we used the SSP grids shown in 
Figure[T9]to constrain the mass of the SC candidate, rejecting all 
SC candidates that lie outside all SSP grids by more than ~ 0.2 
mag. The lower panels of Figure [19] show the mass distributions 
derived from the data-model comparisons described. 

Although we cannot give solid constraints on the present 
mass function (PMF) of SCs, it is instructive to study some fea- 
tures of the magnitude-colour diagrams and PMF derived. First, 
a few very blue [V-I ^ 0.2 mag, o ^ 25 mag] SC candidates 
can only be associated to a metal poor and young component, 
since intermediate [Fe/H] and metal rich models do not provide 
the necessary colour spread to match these sources. Moreover, 
the position in the magnitude-colour diagram of these sources 
implies low masses, M $ 10*Mq. This is also recognised by 
the presence of a peak at ~ IO'^Mq in the PMF at [Fe/H] =- 
1.8, while the PMF derived from other [Fe/H] are nearly flat 
belo w ~ 10^ Mq. The e xpected age-dependent mass complete- 
ness (iGieles et al.ll2007 ^ is probably observed in the panels of 
figure [19] too. In fact, only young SC candidates overlap with 
SSP models with masses below ~ 1.5 • IO'^Mq, while no SC can- 
didate matches with the grid of low-mass f ^ 1 Gyr SSP models, 
no matter what [Fe/H] is adopted. A similar case are the red ob- 
jects at y - / ~ 1.2, which can only be associated to a population 
of old, metal rich SCs. Although the sharp drop of the PMFs at 
masses $ 10"* Mq might be an artifact of the sample incomplete- 
ness, it is worth to note that in all cases a turn over of the PMF 
around ~ 10^ Mq is observed. Finally, a linear fit to the mass dis- 
tribution, within the range of ~ 1O**-^M0H- ~ 1O''-'*M0, provides 
a slope that agrees within Icr uncertaint}(3. 

As an additional test we have obtained masses by using the 
average mass-luminosity (M/Lv) ratios derived from the B - V 
versus V-I colour-colour diagram, then converting V-band mag- 
nitudes to masses. More in details, we have constrained the mean 
M/Ly of the SC candidates by the comparison of SC colours 
with SSP models in a broad range of ages (0.03 < t{Gyr) < 14.0) 
and [Fe/H] from -2.3 to -1-0.3 dex. We set a default 0.25 "colour 
distance" to the SSP model for accepting its M/Ly ratij^. The 
upper panel of Figurel20lshows the colour-colour diagram of SC 
candidates, and the three reference SSP models from SPoT. The 
median M/Ly obtained from the data-to-models comparison was 
then used, together with visual magnitude of the SC candidate, to 
estimate the mass. The main difference with this new approach 
is that no a priori assumption is done on the [Fe/H] content. 
Moreover, in the previous approach, at fixed [Fe/H] the mass and 
age of the SC were virtually present in the comparisons shown 
in Figure [T9l (upper panels), while in this case mass comes from 



' The lower limit to ~ IO^'^Mq has been chosen to avoid the 
masses, i.e. magnitudes, most affected by incompleteness. The B-band 
completeness at masses approximatively lower than ~ 10'*-'Mq drops 
quickly below 90%. 

We made several tests using different ranges (from 0. 1 to 0.5 mag, 
also releasing the matching between colours), without any substantial 
change in the PMF slope. 
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Fig. 19. Upper panels: magnitude-colour diagrams of SC candi- 
dates and SPoT SSP models for the three reference metallici- 
ties ([Fe/H] upper quotes in each panel). Models with ages 0.03, 
0.2, 0.5, 1, 5, 10 and 14 Gyr are connected with a line, older 
ages have redder colours. SSPs models with different masses, 
M,or ~ 5 ■ 10^ 1.5 ■ 10^ 6.0 ■ 10^ 1.2 ■ 10\ 1.0 ■ 10''Mo, 
are shown with increasing size symbols, depending on the mass. 
Models are shifted to the distance modulus of NGC 3370. Lower 
Panels: rough estimates of the Present Mass Functions derived 
from the data to model comparison shown in the upper panels. 
A linear fit to histograms in the range 4.5 < log(M/M0) < 6.5 is 
shown with a long-dashed line. The slope of the linear fit is upper 
quoted in each panel. The shaded region marks, approximatively, 
the mass interval corresponding to a photometric completeness 
below -90%. 



the statistical analysis of the M/Ly of SSP models within large 
age and [Fe/H] ranges. The PMF derived using M/Ly ratios is 
shown in Figure|20] Our main interest, here, is to emphasise that 
the two approaches used to derive the PMF provide similar re- 
sults concerning the power law profile. 

The PMF in spiral and irregular galaxies is generally well fit 
by a power-law, although the fit is limited to the young and most 
massive component of clu sters. The typical value for power law 
exponent is a ~ -2.0 (e.g. lZhang & Fal]||1999HBik et al.ll2003h . 
Similarly, the high mass tail of Galactic OCs has a slope a = 
-2.04 + 0.11, with no selection on clus ters age, thus includin g 
objects from few Myr up to several Gyr dBattinelli et al.lll9^4l) . 

Such numbers agree with the slopes obtained from the 
power-law fit to data (shown in the Figures [T9li20l i. ranging be- 
tween -1.3 + 0.4 and -0.9 + 0.2, and corresponding to a power 
law index a - -2.3 -1.9 i n the form dNsc/dM oc M". 

Recently, iLarsenl (l2009h proposed that the power-law fit to 
the MF of SCs might be due to the limited mass range of obser- 
vational data, and suggested that a functional form of the type of 
a Schechter law - characterised by a power law and an exponen- 
tial decline after some characteristic mass Mc ~ 2.1 x 1O^M0 
- approximates quite well some of the best studied SC systems, 
like in M 51 and LMC. However. ILarsenl also warns that it is 
unlikely that such functional form is universal, quoting the cases 
of some merger or starburst galaxies as examples where signifi- 
cantly higher Mc values are required. As a non-interacting spiral. 
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Fig. 20. Upper panel: SC candidates (dots) compared to SSP 
colour predictions from SPoT for [Fe/H] =-1.8, -0.7, 0.0 (shown 
with solid, dashed and dot-dashed lines, respectively). The mod- 
els age ranges from 0.03 to 14 Gyr. Lower panel: The PMF de- 
rived from the estimate of median M/Ly ratios of SC candidates. 



NGC 3370 should fall in the sample of galaxies with a PMF well 
de scribed by a Sc hechter form. Assuming a Schecter function as 
in lLarsenI (|2009|) . a cut-off mass of the order of Mc ~ 10^ Mq 
would be required for the approximate PMF derived here. Such 
high Mc value has been fit to the young SCs in the Antennae 
dJordan et al.l I2007D . but it is highly unlikely in spirals, being 
probably a result of the sample incompleteness at low masses, 
and the rough mass estimates used here. 

3.5. Globular Cluster candidates 

The colour distributions of NGC 3370 clusters, and the com- 
parisons with GGCs, as well as with the cluster system in other 
galaxies, shown in previous sections, has evidenced the presence 



of a red system of objects with properties remarkably similar 
to old globular clusters. In this section we provide a more de- 
tailed study on such GC candidates. To select a sample of can- 
didate GCs we consider only sources with V - I matching with 
the GC system in the Galaxy, in addition sources in the bright- 
est/bluest part of the disk are rejected to avoid contamination 



> 50", j^v ^ 24.0 mag /arc sec ). With such criteria 



we end up with a list of 35 GC candidates. Figure |2T| shows 
the colour histograms (normalised to the maximum bin ampli- 
tude) for GCs in NGC 3370 and, for comparison, for GGCs 
brighter than My > - 5 (GGC data from iHarrislfl 9961) . The me- 
dian colours are B - V - 0.62 + 0.15 mag (0.68 + 0.11 mag), 
and B - / = 1.60 + 0.27 mag (1.58 + 0.19 mag) for NGC 3370 
(MW) star clusters. The median Rg/f is ~ 5.5 with 3.9 pc stan- 
dard deviation. Having in mind the warning of the small sample 
of GC selected, it is possible to associate the similarity between 
the B-V and B - I colour distributions in both galaxies to sim- 
ilar age/metallicity properties of the two GC systems. However, 
some differences seems to be recognised between the two galax- 
ies. First, the colour distributions in NGC 3370 do not appear so 
sharp as in the case of GGC, though this might be a consequence 
of larger photometric uncertainties. Second, both NGC 3370 dis- 
tributions show a tail of GC at bluer colours with respect to the 
MW - more evident in B - /. If the apparent secondary peak 
is interpreted as a metallicity difference then, using the B - I 
colour-metallicity relation by Harris et al. (2006) derived from 
GGCs, we find that the main peak corresponds to a metallicity 
[Fe/H]~ -1.5 dex (as for the Galaxy), while the [Fe/H] for ob- 
jects in the blue tail is 2.3 dex. However, we must emphasise 

that several doubts have been raised in the last few years against 
the simple interpretation of colour bimodali ties in GC systems as 
metallicity bimodalities. Non-linear metallicity-colour transfor- 
mations, in fact, have been proved to be effective in generating 
multimod al colour distributions from unimodal [Fe/H] distribu- 
tions (see lCantiello & Blakeslee 2007, and references therein). 

One other useful quantity related to the GC system is the 
specific frequency, Sn, i.e. the number of GC normalised to the 
host galaxy luminosity. This number provides a useful quan- 
tity related to the for mation mechanisms of the galaxy (see 
iBrodie & Stra der 2006, and references therein). To properly es- 
timate S ff, however, the total number of GC needs to be used. As 
a consequence of the sample incompleteness we cannot reliably 
constrain the specific frequency, but only provide an estimate to 
the lower limit of 5 at. In particular, we make the follow i ng as- 
sumptions: /) we use the results by, e.g.. iRhode & Zepl (l2003l) 
showing that the surface density of GC in spirals follows a power 
law profile; //) a total apparent corrected magnitude myfi ~ 1 1 .7 
mag for NGC 3370 from HyperledcfHI is adopted, and ///) a the 
slope for the density pr ofile of log pqq ~ - 1 -5 log Rgc is fixed 
from the compilation of lKissler-Patid ( Il997h . for a total galaxy 
magnitude My ~ 21.5 mag. By constraining the intercept of 
the density profile from the number of GCs within 50" H- 80" 
- corrected for photometric completeness - we find a lower limit 
of Spf ~ 0.8, to be c ompared with 5 at > 0.5 for the GC sys- 
tems on other spiral s (iKissler-Patidl 19971; iRhode & Zepll2003l; 
IChandaretaLlllOOl . 



4. Summary 

Deep optical VI and shallower B band observations of the face 
on spiral galaxy NGC 3370, drawn from the ACS/HST archive, 
have been used to study the SC system in the galaxy. The char- 
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Fig. 21. Colour histograms of GC candidates in NGC 3370 (up- 
per panels) and of bright GGCs (lower panels). The histograms 
are normalised to the bin with maximum amplitude. 



acteristics of the imaging data available allowed to analyse the 
photometry of SC candidates, and their spatial properties. The 
final list of SC candidates consists of 277 objects, selected ac- 
cording to maximum colour, magnitude, size and position within 
the galaxy. 

The analysis to the sample of SC candidates leads to the fol- 
lowing conclusions: 

(1) the colour distributions show that the SC population is 
composed by different subpopulations. In particular, the differ- 
ences are related to ages, with the main sub-populations repre- 
sented by /) a red system of likely GC candidates, //) a domi- 
nant peak of sources with properties similar to intermediate age 
clusters (similar to OCs with ages from a few 100 Myrs, to few 
Gyrs), and ///) a blue tail of SCs, possibly with ages below ~100 
Myr. Such complex age-structure of the SC system is somehow 
expected in a Sc-type spiral, actively star-forming (as also testi- 
fies the amount of Cepheids observed in this galaxy); 

(2) the comparison of SC candidates colours - selected with 
no a priori assumption on the range of colours allowed - with 
SSP models predictions within a large range of ages (30 Myr 
to 14 Gyr) and [Fe/H] (-2.3 to -1-0.3 dex) shows a good match- 
ing between models and data. A significant exception to this is 
represented by a sample of ~ 20 SC candidates, which system- 
atically lie outside the region of SSP models, showing an ex- 
cess of flux at longer wavelength (1-band magnitude). Such SC 
candidates do not show remarkable peculiarities in terms of spa- 
tial distributions or morphology. One possible origin of their be- 
haviour is that they are typical intermediate age SCs (300 Myr to 
2 Gyr). At these ages, the stochastic presence of TP-AGB stars 
has a dramatic impact on the integrated light. In particular, TP- 
AGBs domin ate the long-wavelength emission o f the stellar sys- 
tem (see, e.g.. lMarastonl2005l : lMarigo et alj2008i and references 
therein), thus the stochastic presence of TP-AGBs might cause 
the observed p osition of these SCs in the colour-colour diagram 
(lRaimondol2009i) : 



(3) the median R^ff of the system is 4.2 with 3.0 pc standard 
deviation. The R^ff distribution shows a power law drop above 
Reff ~ 3.0 pc. The exponent of the power law, rj = 2.35 + 0.22, 
appears normal with respect to other similar galaxies. Moreover, 
a log-normal fit to the distribution of R^ff seems to show a turn- 
over around ~3 pc, although the completeness of the data does 
not allow to give more robust conclusions; 

(4) the Reff shows non-negligible correlations both with the 
colour of the SC candidate, and with the galactocentric distance 
Rge- If colours differences are associated to age differences, as 
suggested in (1), these results point toward an age evolution 
of Reff , with older SCs having larger R^ff . Furthermore, the 
spatial density of clusters, which appears similar to the surface 
brightness profile of the galaxy, might well be a consequence of 
a common dynamical evolution of stars and star clusters in the 
gravitational potential of the galaxy; 

(5) the LF of SC candidates shows a continuous increase 
down to magnitudes where the SC candidate sample is fairly 
unaffected by incompleteness. Such LF is well described by a 
power law with exponent slightly below ~ -2.0 in all bands. 
This functional form agrees to that of other SC systems in spiral 
or irregular galaxies. A comparison of the maximum SC mag- 
nitude with the SFR per unit area shows that the correlation be- 
tween these quantities is like in other spirals. Moreover, some 
properties of the LFs presented here (steepening of the LF to- 
wards brighter magnitudes, and steeper slopes at longer wave- 
lengths) seem to support the idea that the SCLF is originated by 
a truncated mass function, though sample incompleteness ham- 
pers more detailed conclusions; 

(6) by comparing data to SSP models, we tentatively derived 
a crude estimate of the Present Mass Function of the SC system. 
Since our data suffer for the lack of an age-sensitive indicator, 
the PMF derived must be taken carefully as a best guess with 
available data. With this warning in mind, we find that a power 
law fit to the interval of more massive SCs (i.e. objects likely less 
affected by the incompleteness) provides an a ~ -2.3 H- -1.9, in 
agreement to what can be found in literature for similar galaxies. 
Moreover, the PMF seems to show a turn-over around masses of 
the order of ~ IQ^Mq, though this is possibly due to the sample 
incompleteness at lower masses; 

(7) we analysed the properties of a sample of GC candidates, 
selected on the base of their V - I colour, and their positions 
with respect to the centre of the galaxy. The number of GC can- 
didates selected is 35, these show colour distributions similar 
to GGCs, though less peaked around a median value (proba- 
bly due to the larger photometric uncertainties), and with a rela- 
tively more populated tail towards bluer colours. The peak of the 
colour distribution, converted to metallicity using GGCs data, 
gives [Fe/H] 1.5, similar to GGCs; 

(8) due to the complex issue of completeness in these data, 
we cannot provide a reliable estimate to the specific frequency 
of GCs. However, taking advantage of the known approxima- 
tive behaviour of the surface density profile of GC in spirals, we 
estimate a lower limit Sn > 0.8. 

The high resolution of ACS/WFC allowed us to carry on a 
refined study of the structural and photometric properties of the 
SC system in NGC 3370. These data are essential to provide fur- 
ther constraints to the evolution of stars and stellar systems in 
this galaxy. On the other hand, the lack of age-sensitive indica- 
tors, like shorter-wavelength (~ U-band) photometry, is a serious 
drawback for a detailed study of ages and masses of the SC sys- 
tem in this galaxy. The restoration of ACS and the installation 
of WFC3 on the HST, after SM4, will certainly help to fill the 
present lack of data. 
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